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The non- Abelian geometric piiases of tlie robust degenerate ground states were proposed as phys- 
ically measurable defining properties of topological order in 1990. In this paper we discuss in detail 
such a quantitative characterization of topological order, using generic Abelian fractional quantum 
Hall states as examples. We find a way to determine even the pure U{1) phase of the non- Abelian 
geometric phases. We show that the non-Abelian geometric phases not only contain information 
about the quasi-particle statistics, they also contain information about the chiral central charge Ac 
mod 24 of the edge states, as well as the Hall viscosity. Thus, the non-Abelian geometric phases 
(both the Abelian part and the non-Abelian part) provide a quite complete way to characterize 2D 
topological order. 
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I. INTRODUCTION 

Topological order^'^ is a new kind of order that is be- 
yond Landau symmetry breaking theory, '^'^ and appear 
in fractional quantum Hall (FQH) states. ^'^ It represents 
the pattern of long range entanglement in highly entan- 
gled quantum ground states. The notion of topologi- 
cal order was first introduced (or partially defined) via 
a physically measurable property: the topologically ro- 
bust ground state degeneracies (ze robust against any lo- 
cal perturbations that can break all the symmetries).^ 
The ground state degeneracies depend on the topology 
of the space. Shortly after, a more complete charac- 
terization/definition of topological order was introduced 
through the topologically robust non-Abelian geomet- 
ric phases associated with the degenerate ground states 
as we deform the Hamiltonian.^'^*' It is was conjectured 
that the non-Abelian geometric phases may provide a 
complete characterization of topological orders in 2D,^ 
and, for example, can provide information on the quasi- 
partiicle statistics. Recently, such non- Abelian geometric 
phases were calculated numerically,^^^^'^ and the connec- 
tion to quasi-partiicle statistics was confirmed. 

We like to stress that, now we know many ways to 
characterize topological orders, such as X-matrix,^'^"*"^^ 
F-symbol,^® fixed-point tensor network, etc . But 
they are either not directly measurable {ie correspond to 
a many-to-one characterization), or cannot describe all 
topological orders, or cannot fully characterize topologi- 
cal orders {ie correspond to a one-to-many characteriza- 
tion). So the non-Abelian geometric phases, as a phys- 
ical and potentially complete characterization of topo- 
logical orders, are very important. To understand this 
point, it is useful to compare topological order with su- 
perconducting order. Superconducting order is defined 
by physically measurable defining properties: zero resis- 
tivity and Meissner effect. Similarly topological order 
is defined by physically measurable defining properties: 
topological ground state degeneracy and the non-Abelian 
geometric phases of the degenerate ground states. 
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To gain a better understanding of non-Abelian geomet- 
ric phases of the degenerate ground states, in this paper, 
we are going to study the non-Abelian geometric phases 
of most general Abelian fractional quantum Hall (FQH) 
states in detail. In particular we find a way to determine 
both the non-Abelian part and the Abelian U{1) phase of 
the non-Abelian geometric phases. Through those exam- 
ples, we like to argue that non-Abelian geometric phases 
form a representation of SL{2, Z), which is generated by 
S,T satisfying 5^ = (ST)^ = 1. Such a representation 
contains information on the quasi-particle statistics as 
well as the chrial central charge Ac mod 24. 

It is generally believed that the quasi-particle statis- 
tics and the chiral central charge completely determine 
the topological order. Thus the non-Abelian geometric 
phases provide a quite complete way to determine the 
topological order. In this paper, we will calculate and 
study the non-Abelian geometric phases for a generic 
Abelian FQH state. First, we will give a general intro- 
duction of topological order and FQH states. 




FIG. 1. A particle and its quantum wave on a circle. 



All discovered FQH states have such densities that the 
filling factors are exactly given by some rational num- 
bers, such a.s ly = 1,1/3,2/3,2/5,.... FQH states with 
simple rational filling factors (such as = 1, 1/3, 2/3, ...) 
are more stable and easier to observe. The exactness 
of the filling factor leads to an exact Hall resistance 
Pxy — J^^^-jk- Historically, it was the discovery of ex- 
actly quantized Hall resistance that led to the discovery 
of quantum Hall states. Now the exact Hall resistance 
is used as a standard for resistance. 



II. TOPOLOGICAL ORDERS IN FQH STATES 
A. Introduction to FQH effect 

At high enough temperatures, all matter is in a form 
of gas. However, as temperature is lowered, the motion 
of atoms becomes more and more correlated. Eventually 
the atoms form a very regular pattern and a crystal order 
is developed. With advances of semiconductor technol- 
ogy, physicists learned how to confine electrons on a inter- 
face of two different semiconductors, and hence making 
a two dimensional electron gas (2DEG). In early 1980's 
physicists put a 2DEG under strong magnetic fields (~ 
10 Tesla) and cool it to very low temperatures '^K"), 
hoping to find a crystal formed by electrons. They failed 
to find the electron crystal. But they found something 
much better. They found that the 2DEG forms a new 
kind of state - Fractional Quantum Hall (FQH) state. ^'^^ 
Since the temperatures are low and interaction between 
electrons are strong, the FQH state is a strongly corre- 
lated state. However such a strongly correlated state is 
not a crystal as people originally expected. It turns out 
that the strong quantum fluctuations of electrons due to 
their very light mass prevent the formation of a crystal. 
Thus the FQH state is a quantum liquid. ^'^ 

Soon many amazing properties of FQH liquids were 
discovered. A FQH liquid is more "rigid" than a solid 
(a crystal), in the sense that a FQH liquid cannot be 
compressed. Thus a FQH liquid has a fixed and well- 
defined density. More magic shows up when we measure 
the electron density in terms of filling factor defined as 
a ratio of the electron density n and the applied magnetic 
field B 

_ nhc density of electrons 

eB density of magnetic fiux quanta 



B. Dancing pattern — a correlated quantum 
fluctuation 

Typically one thinks that only crystals have internal 
pattern, and liquids have no internal structure. How- 
ever, knowing FQH liquids with those magical filling fac- 
tors, one cannot help to guess that those FQH liquids 
are not ordinary featureless liquids. They should have 
some internal orders or "patterns" . According to this 
picture, different magical filling factors can be thought 
to come from different internal "patterns". The hypoth- 
esis of internal "patterns" can also help to explain the 
"rigidness" (ie. the incompressibility) of FQH liquids. 
Somehow, a compression of a FQH liquid could break 
its internal "pattern" and cost a finite energy. However, 
the hypothesis of internal "patterns" appears to have one 
difficulty - FQH states are liquids, and how can liquids 
have internal "patterns"? 

Theoretical studies since 1989 did conclude that FQH 
liquids contain highly non-trivial internal orders (or in- 
ternal "patterns"). However these internal orders are 
different from any other known orders and cannot be 
characterized by local order parameters and long range 
correlations. What is really new (and strange) about the 
orders in FQH liquids is that they are not associated with 
any symmetries (or the breaking of symmetries), and 
cannot be characterized by order parameters associated 
with broken symmetries. This new kind of order is called 
topological order, and we need a whole new theory to 
describe the topological orders. 

To gain some intuitive understanding of those internal 
"patterns" (ie topological orders), let us try to visualize 
the quantum motion of electrons in a FQH state. Let us 
first consider a single electron in a magnetic field. Under 
the infiuence of the magnetic field, the electron always 
moves along circles (which are called cyclotron motions) . 
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In quantum physics, only certain discrete cyclotron mo- 
tions are allowed due the wave property of the particle. 
The quantization condition is such that the circular or- 
bit of an allowed cyclotron motion contains an integer 
number of wavelength (see Fig. 1). The above quan- 
tization condition can be expressed in a more pictorial 
way. We may say that the electron dances around the 
circle in steps. The step length is always exactly equal to 
the wavelength. The quantization condition requires that 
the electron always takes an integer number of steps to go 
around the circle. If the electron takes m steps around 
the circle we say that the electron is in the (m -|- 1)*'' 
Landau level, (to = for the first Landau level.) Such a 
motion of an electron give it an angular momentum mh. 
We will call such an angular momentum (in unit of K) the 
orbital spin Sob of the electron. So Sob — fn for electrons 
in the (to -I- 1)*^ Landau level. 

The orbital spin can be probed via its coupling to the 
curvature of the 2D space. As an object with no-zero 
orbital spin moves along a loop C on a curved 2D space, a 
non-zero geometric phase will be induced. The geometric 
phase is given by 



= Sob <P dx'-uji = Sob 



c 



d^xR 



(1) 



Sc 



where Sc is the area enclosed by the loop C. Here R 
is the the Gauss curvature and uji the connection whose 
curl gives rise to the curvature. The relation between 
LUi and R = SijdiUJj is the same as relation between the 
gauge potential and the "magnetic field" b = Sijdiaj. 
On a sphere, we have f , d^xR = An. 

^ ' J sphere 

When we have many electrons to form a 2DEG, elec- 
trons not only do their own cyclotron dance in the first 
Landau level, they also go around each other and ex- 
change places. Those additional motions also subject 
to the quantization condition. For example an electron 
must take an integer steps to go around another electron. 
Since electrons are fermions, exchanging two electrons in- 
duces a minus sign to the wave function. Also exchang- 
ing two electrons is equivalent to move one electron half 
way around the other electron. Thus an electron must 
take half integer steps to go half way around another 
electron. (The half integer steps induce a minus sign 
in the electron wave.) In other words an electron must 
take an odd-integer steps to go around another electron. 
Electrons in a FQH state not only move in a way that 
satisfies the quantization condition, they also try to stay 
away from each other as much as possible, due the strong 
Coulomb repulsion and the Fermi statistics between elec- 
trons. This means that an electron tries to take more 
steps to go around another electron. 

Now we see that the quantum motions of electrons in 
a FQH state are highly organized. All the electrons in 
a FQH state dance collectively following strict dancing 
rules: 



2. an electron always takes odd integer steps to go 
around another electron; 

3. electrons try to stay away from each other (try to 
take more steps to go around other electrons). 

If every electrons follows these strict dancing rules, then 
only one unique global dancing pattern is allowed. Such 
a dancing pattern describes the internal quantum motion 
in the FQH state. It is this global dancing pattern that 
corresponds to the topological order in the FQH state. ^'^ 
Different QH states are distinguished by their different 
dancing patterns (or equivalently, by their different topo- 
logical orders). 

Recently, it was realized that topological orders are 
nothing but patterns of long range entanglements.^ FQH 
states are highly entangled state with long range entan- 
glements. The dancing pattern described above is a sim- 
ple intuitive way to described the patterns of long range 
entanglements in FQH states. 



C. Ideal FQH state and ideal FQH Hamiltonian 

A more precise mathematical description of the quan- 
tum motion of electrons described above is given by the 
famous Laughlin wave function^ 



(2) 



' iyj is the coordinate of the j electron 



where Zj = 

The wave function is antisymmetric when m =odd inte- 
ger and correspond to a wave function of fermions. Such 
a wave function describes a filling factor v = \/m FQH 
state for fermionic electrons. In this paper, we will also 
consider the possibility of bosonic electrons. A v — 1/m 
FQH state for bosonic electrons is described by eqn. (2) 
with TO = even. 

We see that the wave function changes its phase by 
27rTO as we move one electron around another. Thus an 
electron always takes to steps to go around another elec- 
tron in the Laughlin state "^m- We also see that, as 
Zi — )> Zj, the wave function vanishes as {zi — ZjY^ , so that 
the electrons do not like to stay close to each other. Since 
every electron takes the same number of steps around 
another electron, every electron in the Laughlin state is 
equally happy to be away from any other electrons. 

The Laughlin wave function has another nice property: 
it is the exact ground state of the following A^-electron 
interacting Hamiltonian 

^ 1 

1=1 i<j 

where = {xi,yi) is the coordinate of the z**^ electron 
and 



1. all electrons do their own cyclotron dance in the Vm{r) = ^ w;(— )'9'.(5(z)9' 
TO*^ Landau level {ie m step around the circle); .m-i 



iy. (4) 



1 = 1, ■■■ ,m-l 
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To see ^'m({^i}) ^o be the ground state of the above 
Hamiltonian, we note that aU electrons in ^m({-Zi}) ^-^^ 
in the first Landau level. Thus we can treat the kinetic 
energy —J2f=i k^^^i ~ a constant. We also 

note that average total potential energy of 5'm({2i}) van- 
ishes: 

Vtot = J Y[d^z^ *:„({zJ)^K.(z,-Zj)*™({zJ) =0. 

i i<j 

This is because the function d^.^^.'^m still have m — I 
order zero as zi — >• Zj. Thus S{zi — Zj)dl._^.'i'„i = is 
I < m. When vi > 0, the potential energy is positive 
definite. Thus 'ifm{{zi}) is the exact ground state of the 
Hamiltonian (3). 

If we compress the Laughlin state, some m*^ order ze- 
ros between pairs of electrons become lower order zeros; 
so some electrons take less steps around other electrons. 
This represents a break down of dancing pattern (or cre- 
ation of defects of topological orders), and cost finite en- 
ergies. Indeed, for such a compressed wave function, the 
average total potential is greater than zero since the order 
of zero in the compressed wave function can less or equal 
to the I in the potential (4). Thus the Laughlin state ex- 
plains the incompressibility observed for v = l/m FQH 
states. 

Since the ground state wave function of the Hamil- 
tonian Hm (3) is known exactly and describe a FQH 
state, we call Hm an ideal FQH Hamiltonian. Its ex- 
actly ground state \E'm is called ideal FQH states. From 
the ideal FQH wave function 4"^, we can clearly see the 
dancing pattern of electrons. 

D. Pattern of zeros characterization of FQH states 

In the above, we have related the dancing patterns (or 
more generally, patterns of long range entanglements) to 
the patterns of zeros in the ideal FQH wave functions. 
Such an idea was developed further in Ref. 24. (see also 
Ref. 25-30.) 

A set of data called the patterns of zeros^"' was intro- 
duced to characterize the dancing patterns (or pattern of 
long range entanglements). By definition, the pattern of 
zeros is a sequence of integers {50,0 = 2,3,4, ...}, where 
Sa is the lowest order of zeros in the ideal wave function 
as we fuse a electrons together. Naturally a pattern of 
zeros must satisfy certain consistent conditions so that 
it does correspond to an existing FQH wave functions. 
The pattern of zeros that satisfy the consistent condi- 
tions can be divided into two classes: quasi-periodic so- 
lutions with a period n and non quasi-periodic solutions. 
The non quasi-periodic solutions may not correspond to 
any gapped FQH states, so we only consider the quasi- 
periodic solutions. Those solution is said to satisfy the 
n-cluster condition. This quasi-periodicity reduces the 
infinite sequence {Sa,a = 2,3,4,...} into a finite set of 
data {n; m; 5*2, • • • , 5„}. 



It turns out that the pattern of zeros data 
{n; to; S'2, • • ■ ,S'n} can characterize many non-Abelian 
FQH states. ^^''^^ We can calculate the filling fraction, the 
ground state degeneracy, the number of types of quasi- 
particles, the quasi-particle charges, the quasi-particles 
fusion algebra, etc from the data {n; to; S2, • • • , 5„}.'^'^"'^^ 
However, in this paper, we will concentrate on Abelian 
FQH states, which can be characterized by a simpler K- 
matrix. 



E. The {K, q, s) characterization of general Abelian 
FQH states 

Laughlin states with filling fraction v — l/m rep- 
resent the simplest QH states (or simplest topological 
orders). FQH states also appear at filling fractions 
1/ — 2/3,2/5,3/7, .... So what kind of dancing pattern 
(or topological order) give rise to those QH states? 

V = \/m Laughlin states as simplest FQH states con- 
tain only one component of incompressible fluid. One 
way to understand certain more general FQH states 
(called Abelian FQH states) is to assume that those 
states contain several components of incompressible fluid. 
The filling fraction v = 2/3,2/5,3/7,... FQH states are 
all Abelian FQH states and can be understood this way. 

The dancing patterns (or the topological orders) in an 
Abelian FQH state can be described in a similar way by 
the dancing steps if there are several types of electrons 
in the state (such as electrons with different spins or in 
different layers). Each type of electrons form a compo- 
nent of incompressible fluid and the Abelian FQH state 
contain several components of incompressible fluid. It 
turns out that even for an FQH state formed by a single 
kind of electrons, the electrons can spontaneous separate 
into several "different types" and the FQH state can be 
treated as a state formed by several different types of 
electrons. ^^'^^ 

With such a understanding, we find that the dancing 
patterns in an Abelian FQH state can be characterized 
by an integer symmetric matrix K and an integer charge 
vector q. 15. 16,37,38 q^j-^g dimension of K, q, s is k which 
corresponds to the number of incompressible fluids in the 
FQH state. An entry of K, K^'' , is the number of steps 
taken by a particle in the I**^ component to go around 
a particle in the J**^ component. We note that when 
I = J, the order of zero Kn depends on the statistics 
of electrons. For bosonic electrons, Kjj — even and for 
fermionic electrons, Ku = odd. 

An entry of q, qi, is the charge (in unit of e) carried by 
the particles in the /'^ component of the incompressible 
fluid. An entry of s, sj, describe how the /"^ component 
of the incompressible fluid couple to the curvature of the 
2D space. It is roughly given by the orbital spin of the 
particles in the component. Or more precisely, if the 
particle in the 7*^ component are in the to*^ Landau level 
sj ~ m — 1 + ^^^.'^^ The term to — 1 is the orbital spin 
of the particles in the /"^ component. There is also a 
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"quantum" correction that is beyond semi-classical 
picture of electron dancing step. 

All physical properties associated with the topological 
orders can be determined in term of K and q. For exam- 
ple the filling factor is simply given hy v — q^K~^q. 

In the {K, q, s) characterization of QH states, the v = 
1/m Laughlin state is described by K = m, q = I, and 
s — m/2, while the i' = 2/5 Abelian state is described 

hyK^(l 2^,q=Q,ands=(^^/^). 

The topological order described by {K, q) have a pre- 
cise mathematical realization in term of ideal multilayer 
FQH wave function: 

*-={nn(-/--/)""}x 

where /, J — 1, • • • , t label k different layers and zj is 
the coordinate of the z*'' electron in the J**^ layer. For 
such a multilayer system, electrons in each layer naturally 
from a component of incompressible fluid. Since those 
electrons all carry unit charges, so = 1. The exponents 
(or the order of zeros) Kij form the symmetric integer 
matrix K — [Kjj). The topological order in the FQH 
state (5) is described hy K — {Kjj), qj = 1, and s/ = 
Ku/2. 

We would like to point out that there exists a Hamil- 
tonian (of a form similar to that in eqn. (3)) such that 
the wave function 'i>K is its exact ground state. Thus 
is an ideal FQH state. The dancing pattern and its 
{K, q, s) characterization can be seen clearly from the 
ideal FQH wave function ^'j^. 

It is instructive to compare FQH liquids with crys- 
tals. FQH liquids are similar to crystals in the sense 
that they both contain rich internal patterns (or inter- 
nal orders). The main difference is that the patterns in 
the crystals are static related to the positions of atoms, 
while the patterns in FQH liquids are "dynamic" associ- 
ated with the ways that electrons "dance" around each 
other. (Such a "dynamic" dancing pattern corresponds 
to pattern long range entanglements.^) However, many of 
the same questions for crystal orders (and more general 
symmetry breaking orders) can also be asked and should 
be addressed for topological orders. 

F. Chern-Simons effective theory for Abelian FQH 
states 

We know that the low energy effective theory for a 
symmetry breaking order is given by Ginzburg-Landau 
theory. What is the low energy effective theory of topo- 
logical order? One of the most important application of 
the {K, q, s) characterization of Abelian FQH states is 
to use the data {K, q, s) to write down the low energy 
effective theory of the Abelian FQH states. 



There are two types of low energy effective theories for 
Abelian FQH states: Ginzburg-Landau-Chern-Simons 
effective theory'*""'*^ and pure Chern-Simons effective 
theory^'^^''^^. The pure Chern-Simons effective theory 
is more compact and more closely related to the {K, q, s) 
characterization. The detailed derivations of pure Chern- 
Simons effective theory has been given in Ref. 39 and 43 
and I will not repeat them here. We find that an Abelian 
FQH states characterized by (K, q, s) is described by the 
following Chern-Simons effective theory 

C^-T^^ai^d.ajxe^''^ (6) 

where • • • represents higher derivative terms [such as the 
Maxwell term {d^aiv — d^ai^Y]. It is important to note 
that the scales of the gauge fields a/^, are chosen that the 
gauge charges of the gauge fields a/^ are all quantized 
as integers. Thus we cannot redefine the gauge fields 
fl/p ^ X]j Uijaj^ to diagonalize K. 

The Chern-Simons effective theory allows us to cal- 
culation all the topological properties of Abelian FQH 
state. The filling fraction is given hy ly = q^K~^q. The 
quasi-particle excitations are labeled by integer vectors 
I. A component of Z, Ij, corresponds to the integer gauge 
charge of a/^. In fact, a quasi-particle labeled by I is 
described by the Following wave function 

i,I 

The electric charge of the quasi-particle I is given by Q = 
q^ K~^l and the statistics is given by 6* = irl^ K~^l. 

To compare with some later discussions, we introduce 
another way to label the quasi-particles: we use a rational 
vector ct = K~^l to label the quasi-particle labeled by 
I in the previous scheme. A quasi-particle labeled by a. 
carries an electric charge Q = q^a and has a statistics 

We note that, in the new labeling scheme, when a is an 
integral vector a. = n, the corresponding quasi-particle 
excitation carries an integral electric charge Q — q^n. 
The statistics of the excitations also agrees with the 
statistics of q^n electrons. So such an excitation can 
be viewed as a bound state of a integral number of elec- 
trons. We may view quasi-particles that differ by sev- 
eral electrons as equivalent. Thus a. is equivalent to 
ol' = a + n. The equivalent classes of cx correspond 
to different types of fractionalized excitations. They cor- 
respond to a's that lie inside of the k dimensional unit 
cube. Since the number of a's inside of the k dimen- 
sional unit cube is |det(/ir)|, so there are |det(iir)| — 1 
different types of the fractionalized quasi-particles. (Note 
cy. — corresponds to the trivial unfractionalized excita- 
tions which correspond to bound states of electrons.) 
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FIG. 2. (Color online) The plane represents the space of pos- 
sible Hamiltonians. The solid lines are the phase transition 
lines that separate phases with different topological orders. 
The solid dots represent the ideal Hamiltonians whose ground 
state can be characterized by {K,q,s). The dash line repre- 
sents the real Hamiltonians that can be realized in experi- 
ments. 



G. The {K,q,s) characterization is not physical 
enough 

As we have stressed that the key to understand topo- 
logical order is to find a characterization of topological 
order. For the real Hamiltonian of a FQH system 



^ 1 
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(7) 



the ground state wave function is very complicated. We 
can hardly see any pattern from such a complicated wave 
function. 

What we did above is to deform the Hamiltonian into 
an ideal Hamiltonian which have a simple ideal ground 
state wave function. We find that the ideal ground state 
wave function can be characterized by {K,q,s). Then 
we say that all the states in the same region of the phase 
diagram that contains the ideal state are characterized 
by the same {K,q,s) (see Fig. 2). This allows us to 
characterize topological order in the ground state of real 
Hamiltonian. 

Such a understanding also indicates that the {K, q, s) 
characterization of topological order is not very direct. 
We cannot determine K even if we know the ground 
state wave function of a (non-ideal) Hamiltonian exactly. 
In other words, K is not a directly measurable physical 
quantity. We need to rely on the deformation to ideal 
Hamiltonian to have a {K, q, s) characterization of topo- 
logical order. This raises a possibility that two different 
sets, {Ki,qi,si) and {K2,q2,S2), may characterize the 
same topological order, since the two ideal FQH state 
characterized by {Ki,q^, Si) and {K2, q2, •S2) may belong 
to the same phase (see Fig. 2). (ie we can deform one 
ideal FQH state into the other smoothly without encoun- 
tering a phase transition.) Indeed, through a field redef- 
inition a/^ — >■ Wjjajf^ in the Chern-Simons effective 
theory (6), we see that, if {Ki,qi,si) and {K2,q2, S2) 
are related by"^^ 



Wqi = Wsi 



S2, W^KiW = K2, W€SL{k,Z), 

(8) 



then {Ki, q^, Si) and {K2, q2, ^2) describe the same FQH 
states. Note that for W G SL{k, Z) the field redefinition 
fl/p — > J2j Wijaj^ does not change the quantization con- 
dition of the gauge charges {ie all gauge charges remain 
to be integers). 

However, the equivalence condition (8) is not complete. 
It was point out in Ref. 44 that even K matrices with 
different dimensions may describe the same FQH state. "'^ 
Right now, we believe that all Abelian FQH states are 
described by {K,q,s). But {K,q,s) is a many-to-one 
labeling of Abelian topological order. It is not clear if we 
have found all the equivalence conditions of {K,q,s)'s. 

The dancing pattern and the related {K, q, s) descrip- 
tion is an intuitive way for us to visualize topologi- 
cal order (or pattern of long range entanglements) in 
(Abelian) FQH states. But, due to its many-to-one prop- 
erty, {K, q, s) is not something that can be directly mea- 
sured in experiments or numerical calculations. So such a 
description of topological order, although very pictorial, 
is not very physical. A more rigorous and more physical 
(partial) definition of topological order^'^ is given by the 
topology-dependent ground state degeneracy which can 
be measured in numerical calculations. 

In this paper, wen will concentrate on the physical 
characterization of topological orders. We will concen- 
trate on the ground state degeneracy of FQH states on 
compact spaces. We obtain more physical data to char- 
acterized the rich topological orders in FQH states, we 
will also study the non-Abelian geometric phase associ- 
ated with the degenerate ground states as we deform the 
Hamiltonian of the FQH systems. First, to understand 
the ground state degeneracy of FQH states, in the next 
section, we will study FQH state on torus. 



III. FQH STATES ON TORUS 

A. One electron in a uniform magnetic field 

One electron in a uniform magnetic field B is described 
by the following Hamiltonian 



H 



where 



[A,, Ay) = i-By,0). 



(9) 



(10) 



The ground states of eqn. (9) are highly degenerate. 
Those ground states form the first Landau level and have 
the following form 



^'(a;,?/) = «'(z) = /(z)e 



(11) 



where the function f{z) is a holomorphic function of com- 
plex variables z = x + iy 

First let us consider the symmetry of the Hamiltonian 
(9). Since the magnetic field is uniform, we expect the 
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translation symmetry in both x- and y-directions. The 
Hamiltonian (9) does not depend on x, thus 

But the Hamihonian (9) depends on y and there seems 
no translation symmetry in y-direction. However, we do 
have a translation symmetry in y-direction ounce we in- 
clude the gauge transformation. The Hamiltonian (9) is 
invariant under y ^ y + dy transformation followed by a 
QidyBx+ici> [7(i)_gauge transformation: 

dyy ^yy ' y ' 

In general 

TlHTa^H, j,^^ ^idyBx+ifd^dy^d-a_ ^^2) 

where we have chosen the constant phase in Td 

B 



2 dxdy 



(13) 



to simplify the later calculations. The operator Ta is 
called magnetic translation operator. So the Hamiltonian 
(9) does have translation symmetry in any directions. 
But the (magnetic) translations in different directions do 
not commute 

Td.Td, = e'^('^i-''=--'^i-'^^-)Td,rdi. (14) 

and momenta in x- and y-directions cannot be well de- 
fined at the same time. 

We note that Td acts on the total wave function 
f{z)e^2-y . It is convenient to introduce 

B , 2 B ^ 2 



— nTV T.p-2-V 



= e 



— QidyB{z + 



I Q2 y Q 2 



<ix + idy , 



r.d9 



that act on the holomorphic part f{z). 



B. One electron on a torus 

Now, let us try to put the above system on a torus by 
imposing the equivalence or periodic conditions: z ^ z-\-l 
and z ^ z + T where r = r^, -I- ir^ is a complex number. 
Using the magnetic translation operator (12), we can put 
the system on torus by requiring the wave functions to 
satisfy 

7^(i,o)*(a;,y) = *(x,y), r(,^,,^)*(x, y) = ^{x,y); 

(15) 

eqn. (15) gives us 



^'(a; + l,y) = *(a;,y), 

*(a; + T^,y + ry) = *(a;,y)e" 



(16) 



which are called quasi-periodic conditions on the wave- 
function. In order for the above two conditions to be 
consistent, the strength of magnetic field B must satisfy 

Q—iTyB{x + l) _ Q—'lTyBX 



B = 2tt- 



(17) 



for an integer iV^. is the total number of flux quanta 
going through the torus, eqn. (15) implies that 

fimf{z) = /(z), %..,^^)}(z) = /(z); (18) 



or 



/(z + l) = /(z) 

f{z + t)= /(^)c-i^'^^*-i27rJV^^ 



The above can be written as 

/(z-fa-|-6r) /(z)e-'^''^'^*^'2''''^*^ a, 6 = integers. 

(19) 

The holomorphic functions that satisfy the above quasi- 
periodic conditions have a form 



f{z)^Y.^{n)e-^ 



i ^^fL + i 2TTnz 



(20) 

We see that there are linearly independent holomor- 
phic functions that satisfy eqn. (19). So there are 
states in the first Landau level on a torus with iV0 flux 
quanta. Note that one can easily see that /(z+1) = f{z). 
To show f{z + t) = /(^)e-*^'^^*-'2^^*^, we note that 



fiz) = 

n 

n 

Therefore 



wf" +i27rnz^_i2^T-„ 



/(z + r) = e'^'^^'*-'2^^'*("+^)^x(")e 

n 

n 

A particular choice of the linearly independent 
holomorphic functions can be obtained by choosing xl*^) 
in eqn. (20) as 

Xin) ^lifn% N^^l, x{n) = if n % iV^ ^ Z, 

where a%b = a mod b, and I ^ 0, - ■ ■ , N^f, — 1. We obtain 

/('^(z|t) = 0,/^r^(iV^z|iV^T) =2^e * , 

n 

(21) 

All holomorphic functions /(z) that satisfy the condition 
(19) are linear combinations of the above functions. 
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C. Translation symmetry on torus 

The one-electron system on the torus, eqn. (9), has a 
translation symmetry generated by the magnetic trans- 
lation operator Td (12). But to be consistent with the 
quasi periodic condition (15), the allowed translations 
must commute with T(^i,o) and T(^^ t^). Those allowed 
translations are generated by two elementary translation 
operators 



(22) 



The corresponding operators that act on f{z) are given 

by 



^1 - '^(M- 



(r^N-^TyN-^) 



(23) 



We note that 



T1T2 = e'^^/^^TaTi, fifa = e'^'^/^^fifi. (24) 

Such an algebra is called Heisenberg algebra. Since Ti 
and T2 commute with the Hamiltonian (9), the degener- 
ate eigenstates of eqn. (9) form an representation of the 
Heisenberg algebra. The Heisenberg algebra has only 
one irreducible representation which is dimensional. 
Thus the Hamiltonian (9) must have N^f, fold degeneracy 
for each distinct eigenvalue. Those TV^ fold degenerate 
states form a Landau level. In the first Landau level, 
those states are given by eqn. (21) and they satisfy 

A* 4, 

f.fHAr) = /«(. + ^|r)e'^-^'^-^ = /('+i)(z|r) 

(25) 

Note that from eqn. (21), we find 



i^{N4,n+lY+\2TT{N^n+l){z+^) 



= e ' E e' 



/('+i)(z|r)e 



- i ^ — i 27rz 



D. Many electrons on a torus 

A many electron system in a magnetic field is described 
by eqn. (3). Assuming the magnetic field is given by 
eqn. (17) and If the filling fraction v < 1 and the inter- 
action is weak, all the electrons will in the first Landau 
level. The many-electron wave function will have a form 



E2 
Vi 



where we have assumed that the magnetic field is given 
by eqn. (17). Here = a;^ + ij/i is the coordinates for 
the i**^ electron. 



Assuming the space is a torus defined by z '-^ z -I- 1 and 
z ~ z -|- T, the many-electron wave function ^'({zi}) on 
torus must satisfy 

-^{- ■ ■ ,x^ + l,y^,■ ■ ■) ^-^{- ■ ■ ,x^,yi, ■■■) 
^{■- ■ ,Xt + T^.Hi + Tj^, •••)==*(••• ,a;j,yi, ••• )x 

Q-i2-nN^Xi-i-KN^T^ 

(26) 

for each individual i. This leads to a condition on f{{zi\): 

/(•••, z,+a+6r, ...) = /(..., z„ • • .)^-~^r.b-N,-i2.bN,.. ^ 

(27) 

where a,b — integers. 

If the interaction potential between the electrons is give 
by Vm in eqn. (4), than 5'({xi,?/i}) is the exact zero en- 
ergy eigenstate of eqn. (3) if f{{zi}) has a m**^ order zero 
between any pair of electrons: 

/(••• ,z„---) = U.-..0((z.-z,)™). (28) 

The holomorphic function /({z^}) that satisfy both 
eqn. (27) and eqn. (28) exist only when the number of 
electrons N satisfy N < N^/m. When N — N^jj/m, such 
f{{zi]) has a form 

/({zO) = /,(z)/.({zj), /.({z,}) = n [e,M - ^j|t)]" 

(29) 

where Z = ^Zi. The above f({zi\) can be viewed as 
the Laughlin wave function Y\i^j(zi — Zj)"^ on torus [note 
that 011 (zi — Zj\T) ^ Zi — Zj for small Zi — zj]. 
From eqn. (A4) we find 

fr{--- ,Zi + a + bT,---) 
^(_yn(iV-l)(a+6)^^(... ^ . . . ) 

-imTTrb'^{N-l)-i2mTib[J2i{zi-Zj)] 

e ■> , 

fr{' ' ' J Zi, ■ ■ ■) 

^- irub^ i2TvbN^Zi ^_-^Tn{N ~l){a+b) ^imTTrb'^ +'i2rmrbZ 

We see that f{{zi}) will satisfy the quasi periodic condi- 
tion (27), if fc{Z) satisfies 

fciZ + a + br) = /^(^)(_)™(A^-i)('i+b)e-i™^'^'''-i2™''&^. 

(30) 

There are m linearly independent holomorphic functions 
satisfying the above condition (see eqn. (20)). One par- 
ticular choice of those m linearly independent holomor- 
phic functions is given by (see eqn. (A2)) 

jW(^) = cim7rJ + im7r(Z+i)^^ f^rn{Z + ^ + ^) mT^ , 

if m{N - 1) = odd, 
/W(^) = e^{mZ\mT), if m{N - 1) = even, (31) 
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where I = 0, • • • , m 
find that, for ■m{N^ - 



- 1. Note that from eqn. (A2) we 
1) — even, 



f{N) 



t[^'^ and T2" ' generate the following 



/^('' [Z + a + hT)= Oi/ra {m{Z + a + 6t) tot) 

^— \ T-nrnh^ —1271 Zhrn J [l] ^^ -^ 

for m{N^ — 1) = odd, 
/W(Z + a + 6r) 

m{Z + a + br 



Heisenberg algebra 



(34) 



imTT J + im7r(Z+a+;)r+^)/n 



+ + 0) 



m{Z 



2 + 2) 



TOT 



So, if the filling fraction of the many-electron system is 
V = N/N^ = 1/m, the irreducible representation of the 
Heisenberg algebra will be to dimensional. The eigen- 
states of eqn. (3) are all at least TO-fold degenerate, in- 
cluding the ground states. This is another way to show 
the TO-fold ground state degeneracy of the v — 1/to 
laughlin state. 

The TO ground states in eqn. (31) form a particular 
irreducible representation of the Heisenberg algebra (34). 

Note that f^^' acts only on fc{Z). When m{N - 1) = 
even, we have (see eqn. (25) and eqn. (21)): 



im7r(a+6r) — i TTrmh — i27r(Z-j-^ + ^r)h7n 



i-r 



-b i Tirmb 



We see that there are to different ways to put the v = 
1/to Laughlin wave function ni<j(-^j ~ -^i)™ '^^ ^ torus 
and keep the m}^ order zeros. Those m wave functions 
on torus correspond to the to exact ground state of the 
Hamiltonian (3). Thus the = 1/to Laughlin state has 
TO-fold degeneracy on torus. 

We like to stress that the to*^ order zero is a local 
condition (a local dancing rule). From the above dis- 
cussion, we see that such a local condition leads to the 
TO-fold ground state degeneracy on torus which is global 
and topological property. So the discussion in this section 
represents the rigorous mathematical theory that demon- 
strate how local dancing rules (the to"^ order zero) deter- 
mine the global dancing patterns and number of global 
dancing patterns (the ground state degeneracy). In gen- 
eral, if we want to put a holomorphic function with to*** 
order zeros on a genus g Rcimann surface, we will get 
TO^-fold degenerate ground states. 



fW/(0(z|,) = ;(0(^ + 



N 



N 



ff Vi'^(^lr) = e"^^+'^/f (^+ --|r) 
^/i'+i)(Z|r) 
When m{N — 1) = odd, we have: 



(35) 



fr/i'H^|r) = /i')(^ + 



N 



-e— 7.^'^(^|t) 



N 

n1 



(36) 



Note that for ■m{N — 1) = odd, we have (see eqn. (A3)) 

m{Z 



yi/m 



y\2-iTl/m r{l) 



1 

TO 



E. Stability of ground state degeneracy 

For an TV-electron system, the magnetic translation op- 
erators that commute with the many-body Hamiltonian 
(3) are given by 



dyBxi-\-i —dxdy 



(32) 



The corresponding Td operators that act on f{{zi}) are 
given by 



(33) 



Again, on torus, the allowed magnetic translations 



are generated by 



and 



(N) 



m 

^im7r^ + im7r(ZH 



_ gimTT J-|-imTr(Z+-^ + i) p- i T7r/m- i 27r(Z+ i + i r) , 



^(i + l)/m 



77l(Z 



' e 
2 2' 



^im7r(:^)-iTir/m-i2ir(Z+i + iT) ^(' + 1) 
_ g— iT7r/m— i27rZ + 



(Z) 



The eigenvalues of Ti and T2 are pure phases e"''. 
We see that, in the large limit (with N/N^ = 1/to 
fixed), the phases (p are of order 1. Since Ti and T2 are 
translation operators by distances of order \/N^, the mo- 
mentum differences between different degenerate ground 
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states are of order AA: ~ N^. The electron density is of 
order and the separation between electrons is of order 
Ib ^ l/^/N^- We see that Afc > I/Ib- 

When we add impurities 5H to break the translation 
symmetry, 5H will split the m-fold ground state degener- 
acy of eqn. (3) aX v — 1/m. However, each 5H can only 
cause a momentum transfer of order 1/Ib ^ \/-^(t>- Thus 
in the first order perturbation theory, 5H cannot mix the 
degenerate ground state and cannot lift the degeneracy. 
We need a (^iV^)'^ order perturbation to have a mo- 
mentum transfer of order A/c ~ N^. Such a perturbation 
can only cause an energy splitting of order 



-L/i 



where u is a dimensionless coupling constant to represent 
the strength of the impurity Hamiltonian SH , L = 
represents the linear size of the system, and ^ = — In u 
represents a length scale. We see that the ground state 
degeneracy is robust against any random perturbation in 
the thermodynamical limit. The larger the system, the 
better the degeneracy. 

We can also show that, on genus g Reimann surface, 
the ground state degeneracy of the v = 1/m Laughlin 
state is ? Such a ground state degeneracy is also ro- 
bust against any random perturbations. The existence of 
the robust topology-dependent ground state degeneracies 
implies the existence of topological order. 



IV. GENERIC ABELIAN FQH STATES 
A. Generic Abelian FQH wave functions on torus 

We have seen that the topological order in the Laugh- 
lin state n(^i ~ Zj)™e~ 2" SZfi can be probed by putting 
the state on torus, which results in m-fold degenerate 
ground states. A more generic topological order labeled 
by K and q-^ = (1, • • • , 1) is described a multilayer wave 
function (5). Such a topological order can also be probed 
by putting the wave function (5) on torus. 

In terms of the theta function, the multilayer wave 
function (5) becomes (see eqn. (29)) 



*x = /({^/})e-*^'.^[^-l' 
/({z/}) = U{Z'))fr{{zl}) 

fr{{zl}) = [X{X{ef^'{zl~Z-'^\T)]: 



(37) 



{nn^n"(^f-^/k)} 



on a torus. Here 



(38) 



are the center-of-mass coordinates and /({zf }) satisfies 
a quasi-periodic boundary condition 



/(•••, z/ + a + 6r, ••• ) 



iTiTb'^Ns-i2-KbNsz' 



(39) 



where Nij, is the total number of the magnetic flux quanta 
through the torus. From eqn. (A4), we find 

/,(••• ,z/ + a + 6r,...) (40) 

^ (_)-K.Ha+fc)(_)E..i^...^'^(a+6)j^(... 

where is the total number electrons in the layer. 

We know that cx N^. In large N^f, limit, 

the e"''^'^^ term in eqn. (39) must come from the 
e-iEj^f/./TTTb^A"' ^gj.^ -j^ gqjj^ (-4o)_ -phus A^^ and 

must satisfy 



(41) 



Now eqn. (40) becomes 

/,(••• ,z/ + a + for,...) (42) 

-/,(••• ,z„...)e-'""'''^*-'2.''^*"-x 

^_yKu(a+b)(^_-^Y.j KijN-'(a+b)^iKiiTTTb^ + i27Tb^jKijZ-' 

We see that f{{zf}) will satisfy the quasi periodic con- 
dition (39) if fci{Z^}) satisfies 

/,(••• ,Z' +a + bT,---)=f,{--- ,Z',---)x 



(a+b)(K,i-N^) ~iKi,TTTb^^i2TvbJ2j K,jZ' 



(-) 



The above can be written in a vector form: 



f,{Z + a + 6t) (43) 



where 



= (ATll, 1^22, • ■ ■ ,-^kk), 



(44) 



and = (I,-- - ,1). From now on, we will assume 

A^^ = even, (45) 

and the above becomes 

fc{Z + a + bT) 
— f^(^z^){—)'^^'-°'^^^e^^'^^^^^^^^^^^^'^^ (46) 
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To find the center-of-mass wave function that satisfy 
the above quasi-periodic condition, let us consider the 
following functions labeled by a,ri: 

(Z|r) (47) 

where the vector a satisfies 

Ka — integer vector. (48) 
We find that, for integral vectors a and b, 

/(°''')(Z + a + 6T|r) (49) 

^ ^ gi 7r(n^cx.-\-rj)'^ KT{n-\-a-\-r)) 

_ i Ka — i nb^ Krb ^ 



Thus the distinct functions /'^"■''^ (Z|r) are labeled by a 
in the unit k dimensional cube which will be denoted as 
UC. The number of different a in the unit k dimensional 
cube, a S UC, that satisfy eqn. (48) is |det(i4r)|. Thus 
an Abelian FQH labeled by K has \ (lei{K)\-Jold degener- 
ate ground states. 
Note that 



I 



1 



(n+b+cx+riy KT{n+b+a+rj)+i2Tr(n+a+r)y K-{Z-r)) 



where N = J2i is the total number of electrons. So 
the center of mass motion described Z and the rela- 
tive motion described by zf — z'J totally factorize in the 
'j) ^ ground state wave function. 

Also note that the wave function for the center of mass 
motion ^-^^^(2) = /i"^({Z-f})e"5^(I^../!'.')' form an or- 
thogonal basis 



For bosonic electrons, Kn = even and (— )" — 1. 

Thus, the center-of-mass wave function fc{Z) is given by 



/(°)(Z)=77-«(r)/(°^°)(^) 



(50) 



For fermionic electrons, Kjj — odd and (— )'^^(°+^) ^ 1. 
The center-of-mass wave function ,fc{Z) will be 

/(°)(Z)=ry-«(r)/(«'^"'-/2)(Z) (51) 

Note that here, we include an extra non-zero normaliza- 
tion factor 



(52) 



which is a holoniorphic function of r and satisfies 
tj{t + 1) = e'"/i2r/(r), vi-'^/r) = V^v(.t). (53) 

Each center-of-mass wave function fc°'\z) gives rise to 
a degenerate ground state |a; r) (where t describes the 
shape of the torus). The explicit wave wave for the |a; r) 
state is 



(54) 



f^''\{^'})^f^"H{Z'})fr{{4}) 



K.I i,j 



/ i<j 

From the definition (47), we see that 

f{o'+o.,v)(^Z\T) = (Z|t), a = integer vector. 



(55) 



This can be derived bu considering generalized transla- 
tion operators which translate electrons in one or several 
layers. Because \1/^^(Z) with different a are link by 
generalized translation operators which are unitary op- 
erators, this allows us to obtain (55). A discussion of the 
translations that translate all electrons together is given 
in the next section. Since the center of mass motion and 
the relative motion separate, the many-body wave func- 
tions '^xii^i}) are also orthogonal: 



a/3- 



(56) 



We can always choose a normalization so that 

{^k\^k) = ^olP- But in this case, l^*^) will not be a 
holomorphic function of t. In this paper, we will choose 
a "holomorphic normalization" such that l^*^) are holo- 
morphic function of r and satisfy eqn. (56). 



B. Translation symmetry 

The magnetic translation that acts on the holomorphic 
part f{{zl}) is given by 



(57) 



ij 



Again, on torus, the allowed magnetic translations 
are generated by t[^^ = T', 



(N) 



and Tr 



(N) 



T^'^^ and 72^^ generate the Heisenberg 



algebra (34) where N is the total number of electrons on 

all the layers. Again T^^^ acts only on the center-of-mass 
wave function fc{Z). From eqn. (B8) we find, for bosonic 
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TV 



r) 



electrons, 

= e'2--f^°/i°)(-^k), 

and for fermionic electrons, 



TV,:' 



(58) 



(59) 



AT 



N 

n1' 



friz\r) 

= (-)'^^-^e'2'^^''«/i°)(Z|r), 

= (-)'^'Vi"-*-^"^H^|r). 
where we have used 

K-^I - Af/7V^. 
iV^ = (iVi,... ,iV„), 

= A''/ — the total number of electrons. 



(60) 



(61) 



(62) 



Here Nj is the number of /"^ electrons. 

We know that fc°'\z)fr{{z() for ot g [/C corresponds 
to one of the |det(iir)| degenerate ground states, |q:;t). 
Those degenerate ground states form a representation of 
the Heisenberg algebra (34). We find that, for bosonic 
electrons, 



Tri";r) 



(N) 



(63) 



and for fermionic electrons. 



rfV;r) = (-)'^'''e 



— q„i27rq a 



a:T) 



where we have used q = I and assumed -/Va — even. 



(64) 



topological order, the first thing we need is to have a 
complete characterization of topological order; ie to find 
a labeling scheme that can label all distinct topological 
orders. 

One way to find such a complete labeling scheme is 
to use the non-Abelian geometric phase associated with 
the degenerate ground states. This can be done by 
deforming the mass matrix. So let us consider the fol- 
lowing one-electron Hamiltonian on a torus (X^,X^) ^ 
(X^ + 1, X"^) ^ {X^,X-^ + 1) with a general mass matrix: 



lA 



3 h 



where g is the inverse-mass-matrix and 
{A,,A2) = (-2^iV^X^0) 



(65) 



(66) 



gives rise to a uniform magnetic field with flux quanta 
going through the torus. Let z — x + iy = X^ + X'^t 
where t = + iTy is a complex number. We find that 



X^X^ + T^X^ 



y = TyX' 



and 



9x1 = d^, 



dy = — -dxi + —dxi 



Let us choose the inverse-mass-matrix gij such that, in 
terms of {x,y), the Hamiltonian (65) has a simple form 



2 



2 



'dX^ 



iA^ 



Tx d Id 



Ty dX^ Ty dX^ 



lAyf 



d 



1 . d 



^A,f 



rfdX 

Tx , d 
'r^^dX^ 



'dX^ 
*^2)(^-*Ai) 



We find that the inverse-mass-matrix must be 



V. MODULAR TRANSFORMATIONS OF 
GENERIC ABELIAN FQH STATES 



9{r) = 



(67) 



As we have discussed, the ground state degeneracy 
on compact spaces contain a lot of information about 
the topological order. However, those information is not 
enough to completely characterize the topological order; 
ie FQH states with really different topological orders may 
have the same ground state degeneracy on any compact 
spaces. To obtain a complete quantitative theory for 



and [AxiAy) must satisfy 



T^Ay^{l + ''-%)A,-^A, 



Tx . 1 , 

^Ai + -^A2 
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We find that 



Ay = -^2TTN^y. 



Therefore 



(68) 



which is equal to the Hamiltonian (9) (up to the the 
scaling factor Ty). We see that the Hamiltonian (G5) on 
torus (X^X"^) - (Xi + 1,^2) - (Xi,X2 + 1) with a 
r-dependent inverse-mass-matrix (67) is closely related 
to the Hamiltonian (9) on torus z^z + 1^z + t with 
a diagonal mass matrix. 

The ground state of eqn. (9) has a form 
/(z)e^'^^**' Z'^". This allows us to find the ground 
state of eqn. (65): 



==/(Xi+TX2)e'^^*^(^')' 



^iTrN^y^T^/r 



(69) 



Eqn. (69) relates the one-electron wave function on 
torus (Xi,X2) ^ {X^ + 1,^2) ^ {X\X^ + 1), 
^{X^ , X^), with the one-electron wave function on torus 

z^z + l^z + T, f{z)e~'"^'i'y Z"^". Such a result can 
be generalized to the multilayer Abelian FQH state for 
many electrons. From the Abelian FQH state on torus 
zr^z + l-^z + T (37), we find that the Abehan FQH 
state on torus {X\,X^) - {X^ + 1,X^) - {X\X^ + 1) 
is described by (see eqn. (37), eqn. (50), and eqn. (51)) 

$^({X^^X^^};r) 

= raxi^' + ^})e•-^*^^.^[^^'l^ 

ri{4}) = f"i{Z'})M{z',}) (70) 

We see that the r-dependent inverse-mass-matrix (67) 
leads to a r-dependent ground state wave-functions (70) 
for the degenerate ground states. 

The wave functions (70) form a basis of the degenerate 
ground states. As we change the mass matrix or r, we 
obtain a family of basis parametrized by r. The family 
of basis can give rise to non- Abelian geometric phase"*^ 
which contain a lot information on topological order in 
the FQH state. In the following, we will discuss such a 
non-Abelian geometric phase in a general setting. We 
will use 7 to label the degenerate ground states. 

To find the non-Abelian geometric phase, let us first 
define parallel transportation of a basis. Consider a path 
r(s) that deform the inverse-mass-matrix g(ri) to g(r2): 
ri = r(0) and T2 — r(l). Assume that for each inverse- 
mass-matrix (7[r(s)], the many-electron Hamiltonian on 
torus {X\X^) - {X^ + 1,X'^) ~ {X\X'^ + 1) has /c-fold 
degenerate ground states I7; s), 7 = 1, • • • , k, and a finite 
energy gap for excitations above the ground states. We 



can always choose a basis I7; s) for the ground states such 
that the basis for different s satisfy 



h';s\-7-h;s) 
as 



0. 



Such a choice of basis I7; s) defines a parallel transporta- 
tion from the bases for inverse-mass-matrix g{Ti) to that 
for inverse-mass-matrix g{T2) along the path r(s). 

In general, the parallel transportation is path depen- 
dent. If we choose another path r'(s) that connect ri 
and T2, the parallel transportation of the same basis for 
inverse-mass-matrix ^(ri), \j;s = 0) = I7; s = 0)', may 
result in a different basis for inverse- mass- matrix g(T2), 
I7; s = 1) ^ I7; s = 1)'. The different basis are related by 
a unitary transformation. Such a path dependent unitary 
transformation is the non-Abelian geometric phase. '''^ 

However, for the degenerate ground states of a topolog- 
ically ordered state (including a FQH state) , the parallel 
transportation has a special property that, up to a total 
phase, it is path independent (in the thermal dynamical 
limit). The parallel transportations along different paths 
connecting ri and T2 will change a basis for inverse-mass- 
matrix g{Ti) to the same basis for inverse-mass-matrix 
r(r2) up to an overall phase: I7; s = 1) = e"^|7; s — 1)'. 
In particular, if we deform a inverse-mass-matrix through 
a loop into itself {ie r(0) — r(l)), the basis I7; 0) will par- 
allel transport into I7; 1) — e"^|7;0). Thus, non-Abelian 
geometric phases for the degenerate states of a topo- 
logically ordered state are only path-dependent Abelian 
phases e"^ which do not contain much information of 
topological order. 

However, there is a class of special paths which give rise 
to non-trivial non-Abelian geometric phases. First we 
note that the torus {X\X^) - {X^ + 1,X^) ~ {X\X^ + 
1) can be parameterized by another set of coordinates 



X'2 



d -b 

— c a 



X^ 
X^ 



x^ 
x^ 



a b 
c d 



X'2 



(71) 



where a, b,c,d G Z, ad 
rewritten in vector form 



be 



1. The above can be 



X = MX\ M = 



a b 
c d 



eSL{2,Z). (72) 



{X'^,X'^) has the same periodic condition {X'^,X'^) ~ 
(AT'i + 1,A:'2) ~ {X'\X'^ + 1) as that for {X\X^). We 
note that 



dx2 



d -c\ fdx'i 
-b a J \dx'2 



dx'A ^ fa c\ fdxi 
dx'-^J \b dj \dx2 



The inverse-mass-matrix in the {X^,X'^) coordinate, 
f;(r), is changed to 



d —b\ , ^ I d — c 
-c a]3^^A-b a 
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in the coordinate. From eqn. (67), we find 

that 



I _ I y^r' t' 



with 



, 6 + dr 

T = . 



(73) 



(74) 



M = 



The above transformation t — )■ r' is the modular trans- 
formation. We see that if t and r' are related by the 
modular transformation, then two inverse-mass-matrices 
gir) and gir') will actually describe the same system (up 
to a coordinate transformation). (See Fig. 3) 

Let us assume that the path t(s) connects two r's, 
t(0) and t(1), related by a modular transformation 

- ^ig- We wiU de- 

note t(0) = r and t(1) = t'. The parallel transportation 
of the basis |7;t) for inverse-mass-matrix (7[r(0)] gives 
us a basis |7;t') for inverse-mass-matrix g[T(l)]. Since 
T = t(0) and r' = t(1) are related by a modular trans- 
formation, g[r(0)] and 3[t(1)] actually describe the same 
system. The two basis I7; r) and I7; r') are actually two 
basis of same space of the degenerate ground states. Thus 
there is a unitary matrix that relate the two basis 



|7;t') =C/(M)|7;t) 
C/^y(A/) = (7;r|7';r') = (7; T|f7(M)|7'; r). 



(75) 



Such a unitary matrix is the non-Abelian geometric phase 
for the path t(s). 

Let us consider a modular transformation generated 

by M' = ^,',1 : r' ^ r" = ^^±^. We see that 



c' dl 



b' + d'r' 



C'T' 



{aU + bd') + {cb' + dd')T 
{aa' + be') + [ca' + dd)T 



(76) 



and the transformation r — !• t" is generated by MM' . 
From |7;t") = U{M')\r,T') = C/(M');7(Af)|7; r), we 
find that 

U,y,{MM')^{r,T\j";r") (77) 
7' 

Except for its overall phase (which is path dependent), 
the unitary matrix f/ is a function of the modular trans- 
formation M. From eqn. (77), we see that the unitary 
matrix U form a projective representation of the mod- 
ular transformation. The projective representation of 
the modular transformation contains a lot of information 
(may even all information) of the underlying topological 
order. 

Let us examine (7';T'|7;r) in eqn. (75) more care- 
fully. Let <i>^[{X[}|r] be ground state wave functions 



for inverse-mass-matrix (/[r], and <I>-y[{Xi}|r'] be ground 
state wave functions for inverse-mass-matrix g[T']. Here 
Xi — {Xj' ,X^' ) are the coordinates of the i"^ elec- 
tron in the 7**^ component. Since r and r' are related by 
a modular transformation, $^[{X[}|r] and $-y[{X[}|r'] 
are ground state wave function of the same system. 
However, we cannot directly compare $^[{X[}|r] and 
}|t'] and calculate the inner product between the 
two wave functions as 



ij 



The wave function $.^[{X-}|t'] for inverse-mass-matrix 
g[T'] can be viewed as the ground state wave function for 
inverse-mass-matrix ^[r] only after a coordinate transfor- 
mation. Let us rename X to X' and rewrite $^[{X[}|t'] 
as $^({X^'^}|t'). Since the coordinate transformation 
(72) change r to t' , we see that we should really compare 
$,({X:-^}|r') = %{{M~^Xi}\T') with <i>^({X,[}|r). 
But even $^({Xf}|r) and <^^{{M-'^XI}\t') cannot be 
directly compared. This is because the coordinate trans- 
formation (71) changes the gauge potential (66) to an- 
other gauge equivalent form. We need to perform a C/(l) 
gauge transformation Ug{M) to transform the changed 
gauge potential back to its original form eqn. (66). So 
only $^[{X[}|r] and UG^-y[{M-^Xl}\T'] can be di- 
rectly compared. Therefore, we have 



U^yiM) 



(78) 



= / lld^Xi<i>;i{Xl}\T)UGiM)^yi{M-'Xl}\T') 

which is eqn. (75) in wave function form. Note that t' = 

M12+M22T 
M11+M21T • 

Let us calculate the gauge transformation Ug {M) . We 
note that ^^{{X['^}\t') is the ground state of 

k,I ij = l,2 "^k "^k 

where k — I,-- - ,N labels the different electrons. In 
terms of X' (see eqn. (71)), H' has a form 



H' 



where 



kj i,j=l,2 



d_ 
dXi 



A,\_f d -c\ (A', 
A2 



b a I \A', ' 



a c\ f Ai 



b d \A 



Since {A'^,A'^) = {-2tiN^X''^ We find that 

(ii,i2) = (-27riV0X'2d, 2^N^X'%) 
= [-27: N^{-cX^ + aX^)d, 2ttN^{-cX^ + aX^)b) 
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-2 -1.5 -1 -0-5 O.S 1 1.5 2 



FIG. 3. (Color online) Modular transformations generate by r — >■ r + 1 and r — > — 1/r map r's in one domain to r's in another 
domain, r's in one domain label distinct systems, r's in different domains related by a modular transformation describe the 
same system (up to a coordinate transformation). Each domain has a topology of a sphere minus a point. (The figure comes 
from Wikipedia.) 



Ug{M) will change H' to H: the wave function ^^{{XI}\t): 

I /q . . ^UGiM)^,i{M-'xl}\r') 

= -E2 E (^-i^O^^n^^^'^O ^^yi{Xi}\T)Uy,iM). (80) 

k,I ij = 1.2 '-'^k 

where r' = 4fi4#^. Let 

with (^1,^2) = {-2TrN^X^,0). We find that 
Ug{M;{X,}) ^l[uG{M;Xi), 

k.I 

uc{M;X) = c'^'"'!"-*'*'-*'*"'"-*"-*"*'!, (79) 



Since $^({X[}|r) has a form 
$,({X[}|r) = /,({X/'VrX,'''}|r)e'--^^^.(^^')^ 



Note that at X, = 0, the fixed point of coordinate trans- ^here /^({zflk) is a holomorphic function, we can ex- 
formation = M-^X,, UGiM\X)\x=o = 1. press the action of U{Mt) and C/(Ms) on $^({X[}|t) 



Eqn. (78) can also be rewritten as a transformation 



on as 



U{MT)'S>,{{Xi}\r) = UG{MT)<f,{{Xy Xf•^Xf■^}|r + 1) 

- UG{MT)fj{{{xl'' - xf^') + (t + i)xI'}\t + i)e-(-+i)^* 

and 

UiMs)'f,{{Xl}\r) = C/G(Ms)$,({X^^ ~Xy}\ — ) = C/G(MT)/,({Xf'^ - (Zi)x^^}| ^)e-(^)^* 

T T r 

= e-2-^*5:.,.x^^x^^ f ({^^/^}|Zi)e-i-^.E.,.(^^^)Vr ^ , (|^//^||Zi)e--^*i:.,.(-^)V-e--^^ 

T T 

I 

where z/ = X}'^ + rXf . We see that the action of operator [/(A/) can be expressed as the action of another 
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operator U{M) on /^({X^}|r), where 
We find 



C/(MT)A({z/}|r) = ^({z/}|T + l) 

= fY{{z!}\r)Uy^iMT), 



U{Ms)M{74}\t) ^ /^({z//r}|— )e-'-^*5:.,.K0V^ 



= fy{{z',}\T)Uy^iMs). 



(81) 



We may use this result to calculate U^y (M) . 

In the above, we have assumed that $7(t) form an 
orthonormal basis ($-^(t)|$-^/ (r)) 



In this case, 



U-y-yi (M) is an unitary matrix. However, it is more 
convenient to choose "chiral normalization" such that 
{^^{t)\^^'{t)) cx S-yj' and that ^^(t) is a holomorphic 
function of r. We will choose such a "chiral normaliza- 
tion" in this paper. 

To summarize, there are two kinds of deformation 
loops t(s). If t(0) = t(1), the deformation loop is con- 
tactable [ie we can deform the loop to a point, or in other 
words we can continuously deform the function r(s) to a 
constant function r(s) — r(0) = t(1)]. For a contactable 
loop, the associated non-Abelian geometric phase is ac- 
tually a C/(l) phase U^^' = e^'^S^y. where ip is path 
dependent. If t(0) and t(1) are related by a modular 
transformation, the deformation loop is non-contactable. 
Then the associated non-Abelian geometric phase is non- 
trivial. If two non-contactable loops can be deformed into 
each other continuously, then the two loops only differ by 
a contactable loop. The associated non-Abelian geomet- 
ric phases will only differ by a an overall U{1) phase. 
Thus, up to an overall U{1) phase, the non-Abelian geo- 
metric phases Uj^' of a topologically ordered state are de- 
termined by the modular transformation t ^ t' — ^±*^ . 
We also show that we can use the parallel transportation 
to defined a system of basis ^.y{{Xi}\T) for all inverse- 
mass-matrices labeled by r. 



the non-Abelian geometric phases are determined by the 
center-of-mass wave function only, because the relative- 
motion wave function does not depend on a. If we only 
look at the center-of-mass wave functions <I>^ ^{Z\t), one 
can show that they form a system of orthogonal basis 
(parametrized by r) that are related by the parallel trans- 
formation (up to an overall phase). 

So if we choose the basis for the degenerate ground 
states as |a;T) — $^({X[}|t) (see eqn. (50), eqn. (51), 
and eqn. (54)), then we can use eqn. (81) to calculate the 
non-Abelian geometric phases associated with the two 
generators T : r — > r + 1 and : r — > l/r of the modular 
transformations. 

Let us consider the bosonic electrons. From eqn. (54) 
and eqn. (50), we see that 

/<°n{^/}k) = v-''{r)f^''''\{Z'}\r)x (83) 
Using eqn. (BIO), eqn. (A5), and eqn. (53), we obtain 

(84) 

where TV is a column vector with component Nj and k 
is a column vector with component Ku (see eqn. (62)). 
Therefore, from the definition eqn. (70), we find that 

<i>^{{Xl^',Xf'};T+l) = e''^^*^'.^^^''')^ (85) 
where 



A. Calculation for Abelian states 

For the Abelian FQH states, we have already intro- 
duced a system of basis for the degenerate ground states 
for each inverse-mass-matrix g{T): \a;T) — "I>^({X[}|t) 
(see eqn. (50), eqn. (51), and eqn. (54)). It appears that 
such a system of basis happen to be a system of orthogo- 
nal basis that are related by the parallel transformation 
(up to an overall phase), ie the basis satisfy 



xy = xy + xr, xf^^^xf-^. 



(86) 



(87) 



Using eqn. (B13), eqn. (A6), and eqn. (53), we obtain 
that 



(/3; 



(82) 



One way to see this result is to note that the wave 
function <I>^({JC[}|t) factorizes as a product of center- 
of-mass wave function and relative-motion wave func- 
tion (see eqn. (54)). So up to the over aU U{1) phase. 



/^"H{^/M|-l/r) = (^ 



(N^ KN-k/ N)/2 



X 



^ittZ'^KZ/t 



E 



peuc 



x/ditpO 



'-f^H{4}\r) 
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Note that 



I<J;i,j 
1 

" 2 



I;i<j 



I ,J;i^j I ,i 



(89) 



where we have used iV^ = KjjNj. Therefore, we have 

/(°)(KVr}|-l/r)e-'?^*^-(^^^)^ 

pt^c Vdet(/^) 
From the definition eqn. (70), we find that 
$^({l,^^^lf'};-l/r) = X (90) 



(-) 
where 



„- 121^13 Key. 



E 

/3e(7C 



1 



+TX/-o/r (91) 



(92) 



So for bosonic electrons, the non-Abehan geometric 
phases are described by the fohowing generators of (pro- 
jective) representation of modular group 

U^p{Mt) = eiJ(^"^^-'^"^)e''^(°"^°-A«)j„^ 
Uo.,iMs) = (-)(^-A-iv-.^iV)/2 ejgg- (33) 

Now we consider the fermionic electrons. From 
eqn. (54) and eqn. (51), we see that 

.f^"H{4}\r) = r;-''(r)/(°'^"'-/2)({Z^}|T)x (94) 
/<•/ «,i 



I i<] 

Using eqn. (Bll) and eqn. (A5), we obtain that 
/("H{zf}|T + l) = eA'-(^^^^--^^-'')x 

Using eqn. (B15) and eqn. (A6), we obtain 

/(°)({z//t}| - l/r)e~'^^*^^.-("-)' 

E f'\{zl}\r) 



So for fermionic electrons, the non-Abelian geometric 
phases are described by the following generators of (pro- 
jective) representation of modular group 



e2 



(96) 



VI. SUMMARY 

Let us summarize the results obtained in this paper. 
For an Abelian FQH state (if, q) on a torus, the ground 
state degeneracy is given by |det(i4r)|. The degenerate 
ground states jo;) can be labeled by the k dimensional 
vectors a in the n dimensional unit cube, where k is the 
dimension of K and a satisfy the quantization condition 

Ka. — integer vectors. 

a's that differ by an integer vector correspond to the 
same state. 

Assume that the torus is a square of size 1: (Xi, ^2) ^ 
[Xi + 1,^2) {Xi,X2 + 1) with magnetic flux 
quanta. To simplify our result, we will assume = even 

(see eqn. (45)). The Hamiltonian commutes with t[^^ 

and T2^\ where t}^-* generates translation {Xi,X2) ^■ 

(Ari + -^, ^2) and t!^^^ generates translation (Xi, ^2) 

{Xi,X2 + ^). t[^'^ and t!^^ satisfy the Heisenberg al- 
gebra 



rp{N)rj.(N) 



where N is the total number of electrons. The degenerate 
ground states \a) form a representation of the Heisenberg 
algebra: for bosonic electrons. 



Trie) 



ei2nq «. N 



TfV> = l«- 



and for fermionic electrons, 



K-'q) 



(97) 



Tn«)=(-) 



— /NK q„i2irq a 



a) 



(98) 



For generic inverse-mass-matrix g in eqn. (67), the de- 
generate ground states have a r dependence |a;T). By 
deforming t to t' — . we can obtain a non-Abelian 
geometric phase U{t ^ r'). U{t r')'s form a projec- 
tive representation of the modular group. 

For two generators of the modular group S = U{t ^ 
-1/r) = U{Ms) and T = U{t ^ t + 1) = U{Mt), we 
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find that, for bosonic electrons (see eqn. (93)) 



and for fermionic electrons 



Sep = (-)(^"^^-"^)/2 j2 



v/d^t(Z) 



(99) 



peuc 

and for fermionic electrons (see eqn. (96)) 

Sep = (_)(^"^^-"^)/2 



ap 



-i27T{a^Kp+p'^K.) 



(100) 



where = {K^^ , K'^'^ , ■ ■ ■ . K'^'^) is two times the spin 
vector of the Abelian FQH state. We see that the S- 
matrix contains information about the mutual statistics 
2tioF K(3 between the quasi-particles labeled by a and 
/3, while the T-matrix contains information about the 
statistics ttoFKo. of the quasi-particle labeled by a.}'^^ 
We have calculated the above unitary matrices S", T 
carefully such that the over all [/(I) phase is meanful. We 
have achieve this by choosing the degenerate ground state 
wave functions on torus to be (1) holomorphic functions 
of r and (2) modular forms of level zero (the r — > — 1/t 
modular transformation do not have the \/— Tr factor). 
(1) is achieved by choosing proper gauge for the back- 
ground magnetic field, and (2) is achieved by including 
proper powers of Dedekind rjij) function in the ground 
state wave functions (see eqn. (50) and eqn. (54)). So 
the amibguity of the ground state wave function is only 
a constant factor that is independent of t. Such a con- 
stant factor is a unitary matrix W that is independent of 
T. So the ambiguity of T, S is 



T^WTW\ S~^WSW^. 



(101) 



We note that Ms generates a 90° rotation (see 
eqn. (72)). Also when r = i, the torus has a 90° sym- 
metry. So 5' = i?9oo is the generator of the 90° rota- 
tion among the degenerate ground states when t — i. 



We also note that Ms Mr 
-1 



-1 

1 1 



and {MsMrf 



Q ^, which is a 180° rotation (see eqn. (72)). So 

MsMt generates a 60° rotation. When r = e''^^'^, 
the torus has a 60° symmetry. So U{MsMt) = 
U{Ms)U{Mt) ^ ST Reoo is the generator of the 
60° rotation among the degenerate ground states when 
T = e"^/'^. Wc find that, for bosonic electrons 



iJfiOo = eA-(^"^^— "J^)e-Ai- 



i7^(/3-2a)^i^/3 

J- 1 TTK ^ 



ydet(Z) 



(102) 



^ iTr(Ar^KJV-K^Ar) p-x2 iTTK i iirft^A'^V 



i2Tr (a"^ Kp+p"^ k) 



e2 X 

tt{P+^K-^k)'^ K(P+^K-^k) 



^\n(P-2a-K-^K.)'^ Kp 



(103) 



We have R%qo = i?|oo — 1. 

In S, T, we have a term that depend on KN — 
K^N. Such a term comes from the local relative motion 
(the fr part in eqn. (37)) and is a pure U{1) phase. Other 
parts of T, 5* come from the global center-of-mass motion 
(the fc part in eqn. (37)) and are robust against any 
local perturbations of the Hamiltonian. Since Rqqo = 
Rqqo = 1 and the KN — k^N term is only a pure 
U{1) phase, we see that such a [/(I) phase is quantized 
and is also robust against any local perturbations of the 
Hamiltonian. Therefore, S and T are robust against any 
local perturbations of the Hamiltonian. 

We note that when Nj — mod 24, the U{1) phase 
from N'^KN - k^N is zero (mod 27r). In the following, 
wc will assume Nj = mod 24 and consider only the 
bosonic case. In this case, S, T have simple forms where 
T is diagional and S contains a row and column of the 
same index which contain only positive elements and the 
element at the intersection of the row and the column is 
1. (For more details see Ref. 11-13, and 47.) We can 
make such a row and column to be the first row and the 
first column. In this case, Tn — Ta=a.a=o = c^^^'^'^^^^ 
where Ac = k is the chiral central charge. (The chiral 
central charge is c — c where c is the central charge for the 
right moving edge excitations and c is the central charge 
for the left moving edge excitations.) We see that the 
non- Abelian geometric phases contain information of the 
chiral central charge mod 24. 

We like to remark that if we move r around a samll co- 
tractible loop, the induced non- Abelian geometric phases 
are only a pure U{1) phase. Such a pure C/(l) phase is 
related to the orbital spin^^ or the Hall viscosity. '^^"^'^ 

Although {K, q, s) are not directly measurable, the 

ground state degeneracy and (T;^^^\ T2^\ S", T) are di- 
rectly measurable (at least in numerical calculations). 
Those quantities give us a physical characterization of 
2D topological orders, including the chiral central charge 
mod 24. 
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Appendix B: Multi-dimension theta functions 



Appendix A: Theta functions 

The theta function in eqn. (21) is defined as The Riemann theta function is given by 



9a(z I r) — exp{i7rT(n + a)^ + i27r(r7, + a)z}. 



9a{z I r) has one zero on the torus defined by 
and z ^ z + T. (?q,(z|t) satisfies 



z - z + 1 



,(Z + a + 6T|r) = Qi^r(n+<^f+i2^{n+<^)(z+a+br) 
n 

(n+Q+b)^ + i27r(n+Q)2 

n 

= e^2''""6l„(z|T)e-'^'''''-'2'^^^ a, 6 = integers. 

(Al) 

In particular, we find 

9i,^{m{z + a + 6T)|mT) - 0,/„(™z|mr)e-'"--^'-'2..b™^ 

(A2) 

and 

0,/™(m(z+ -)|mT) = 0i/„,(mz|mT)e'2-'/™^ 
m 

0i/„(m(z + -)|mT) = 0(i+i)/„(mz|mT)e--'^/'"-'2-. 
m 

(A3) 

Let us introduce four variants of the theta function: 
eoMr) = Oo{z\T), 

e,o{z\T) = c^^'^+^'^edz + \T\T), 
eMT) = eo{z + ]^\T), 

e,Mr) = ei'^-+-'(^+5)^o(^ +\ + \r\T). 
^00 7 ^01? E^nd are even and On is odd in z: 

0oo(2|r) = eoo(-^|r), 0oi(^k) - 0oi(-^|t), 
Ow{z\t) = 0io(-z|t), 0ii(z|r) = ^e^i{-z\T). 

From eqn. (Al), we find that 

eii{z + a + bT\T) (A4) 
= 6'ii(z|r)e'^'^''+''^)e~'^''''^^'^''^^+5+i^)''^ 

= (-)''+''6'll(z|T)c-'^'^'''-'2'^^^ a, 6 = integers. 
We also have 

(A5) 



(Bl) 



for any symmetric complex matrix J7 whose imaginary 
part is positive definite. It has the following properties:^^ 



Q{z + a + nb\n) = e 



Q(z\n), (B2) 



for any integer vectors a and b. 
We also have 

e(^|r! + J) = e(^ + ^'^|^|r!), (b3) 

for any symmetric integer matrix J, and 

Q{n-^z\ - n-^) = ^/dei{-\n)e'''^^^''''Q{z\n). (B4) 
These two relation lead to'"'^ 

e([(Cf7 + DfY^z\{An + B)(OT + D)-^) 
= C7\/det(Cf} + i^)e'""^(^^^+^)"'^"e(^|r!), (B5) 
where 

^= (c i^) e^P(2A^,Z), 

diag(A"^C), diag(_B"^£') = even vectors. (B6) 

C-y is an 7 dependent phase satisfying C,^ = 1. 
Let 

FK(Z|r) = Q{KZ\Kt). 

We find that 

Fi^(Z + a + brlr) = e(X(Z + a + br)|i^T) 

for iCa, b = integer vectors. 



and 



Fk{Z\-\It)=Q{KZ\-KIt) 



= Vdet(-iX-iT)e''^^-^ ^^e(TZ|ii:-V) 



v/det(-iX 



v/det(-iX-ir)e''^^-^ ^^x 



E E 



riixin+OL)^ 7< (n+a)+ i 27r(n+Q;) /fZr 



yii(z|T + l) = e3''^6 



qGC/C TiSZ" 

= v/det(-iX-iT)e''^^-^''-^^x 

Qri7ro.^Ko.+ i27ra.^KZrp^(^^^2: + a)\T) 



9lliz\T) [AO) aeuc 

9n(-|--) = -(-ir)i/^e^'^^0n(.|r) (A6) = E Vdet(-i/^-ir)e^'-(^+°)^^'W«)^^^,(,(z + a)|r) 

r r — 



aeuc 
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Let us define (Z|r) as 

/(°''')(Z|r) 

_ gi7r(n+a+r7)^i<'T(Ti+a+77)+i27r(Ti+a+r7)^if (Z-ry) 

neZ" 

_ gi7r(a+rj)^A'r(a+r7)+i27r(a+77)^K(Z-77) ^ 

^ ^ ^iirn^ KTn+i27TTi^ K[Z+((x+ri)r—r]] 

_ giTr(a+77)^Kr(a+r7)+i27r(a+77)^A'(Z-r() ^ 

e(if [Z + (a + J7)r - f]] \Kt). 

_ ^iTT(a+T,)'^ KT{cx+ri)+i2Tr{a+rj)'^ K(Z-ri) ^ 

FK{Z + {a + r])T-r]\T). 

It has the following properties: For any integer vectors a 
and b, we have 

/("■^)(Z + a + 6r|T), 



or 



/(°''')(Z + /^-ia|r) 



(B8) 



e ' e 



and 



/("■'')(Z|t + 1) 
+»7)^is:(T+ 

Fk(Z + (a + 77)(r + 1) - tjIt + 1) 



giir(a+77)^_fs:(T+l)(a+77)+i27r(a+77)^A:(Z-r;) ^ 



+77)-' /i'(T+l)(a+77)+i27r(a+77)-' K(Z-'q) 



Fk{Z + {a + 'q)T + a. - 



_ gi7r(a+77)-' _R'(T+l)(a+77)+ i 27r(a+77) K{Z~ri) ^ 

i^if(Z+(a + J7)r+^^^^|T). (B9) 
where we have used Fk{Z\t) — Fk{Z + cx.\t). Here 

K = {Kii,K22,--)- 

If 77 = and k/2 is an integer vector, we have 

f^"-^°\Z\T + 1) = Q^^<^^K{T+l)o.+ i2^o.^KZp^^^2 ^ ^^1^^ 

= e'''°''-^°/(°^°)(Z|T). (BIO) 
If = we have 



/("■'')(Z|t + 1) 

7r(a+77)^/i'(T+] 

Fk(Z+ (a + r7)r-r7|r), 



gi7r(a+77)-' A (r+l) (0+77)+ i 27r(a+77) ^ ^(^-»?) ^ 



^ ^iTT(Ka+ln)'^ K-\Ka+^K.) j(a,K-^K./2) ^gj^^^j 



We also have 



/("■^)(Z|-1/t) 

i^^K(Z-(a + J7)T-i-J7|-l/T) 

g-T-i i7r(a+?7)^A'(a+r7)+i27r(a+r))'^A'(Z->)) ^ 

^ v/det(-ii\:-iT)x 
/3ec/c 

gri7r(Z-(a+Tj)T-i-r7+/3)^A'(Z-(a+rj)r"i-r7+/3) ^ 

Fk(t(Z- (a + J7)T-i -T7 + /3)|t)x 
^ v/dct(-iA'-ir)x 

gri7r(Z+/3-r7)^A(Z+/3-r7)-i27r/3^A(a+J7)^ 

Fk{tZ - {a + ri) - T-q + r/3|r) 

^ Vdet(-iX-iT)x 
peuc 



gri7r(Z+/3-77)^ A(Z+/3-rj)- i27r/3^ A(a+»7) ^ 

F/aTZ + r(/3-r7)-J7|T). 

where we have used Fk{Z\t) — Fk{Z + q:|t). 
When 77 = 0, we find 



(B12) 



Act, 



/(°^°)(Z|-l/r) 
= ^ Vdet(-i/^-ir)e"-(^+^)"^(^+^)-'2-/3' — X 

/3e(7C 

FK(rZ + r/3|r) 

= Vdct(-iX-lT)c^'^^"^^ ^ e-'2-/3"^oj(/3,0)(^_^|^)^ 



/(«'°)(Z/T|-l/r) (B13) 

_-i27r/3^Aa 

(_i,)./2^i.z-AZ/. ^ e fi^-)i^z\r). 



Note that 



^-i2TTa'^ K-f p+i27r/3^A'7 



^^^^ yd^ yd^ 

Eqn. (Bf3) can also be rewritten as 

/("■°)(Z|t) = (-ir)-'"/2g-i7rZ^AZ/r^ ^g^^^ 
„i2ir/3^Aa 

/atl^c Vdet(A) 
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Whenry = K-'^k/2, we note that Fk{tZ + T{f3-'q)- 
n\T) = Fk{tZ + t{(3 + ^\r)- Thus 

/(«''')(Z|-l/r) 
fseuc 

Fk{tZ + T{p-ri)+r)\T) 
= v/det(-iif-iT)e""^^^^'^x 



Appendix C: Algebra of modular transformations 
and translations 



The translation Td and modular transformation U{M) 
all act within the space of degenerate ground states. 
There is an algebraic relation between those operators. 
From eqn. (80), we see that 

U{M)Td%{{X,}\T) = U{M)%{{X, + d}\T) 
= UG<^>y{{M-'X, + d}\T') 
= C/G$^({Af-i(X, + Md)}|r'). 



Therefore 



or 



e''U{M)Td - TMdU{M). 



(CI) 



/(°-^'''^/2)(Z/t| - 1/r) 



^ (■_j^)K/2gi7rZ^A:Z/Tgi7rK^_ft:-iK/2^ 

E ^ /(/^-^---/^)(Z|r), (B15) 

nt^c Vdet(i^) 



Let us determine the possible phase factor e'^ for some 
special cases. Consider the modular transformation t 
t' = T + 1 generated by 



1 1 



1 -1 
1 



which can also be rewritten as 



/(°'^"'-/2)(Z|t) 



(B16) 



(-ir)-''/^e 



k/2 -i7rZ^KZ/T -i7rK^if-iK./2^ 



„i27r(a^ K.) 

E /(/^^^--/-)(Z/r|-l/r). 
^ti/c Vdct(i^) 



From eqn. (B14), we find that 
/(«'°)(TZ|r) 



(B17) 



/3ec/c \/det(if) 

,i27r/3^i<'Q; 



We first calculate 

U{Mt)T^^^{{X,}\t) = U{Mt)^^{{X, + di}|r) 
= e-'"^*2:.(x?)^$^({M^iX,; + di}|r'). 

where di = (-^,0). We note that 

Thus we next calculate 

TiU{Mt)%{{X,}\t) 
^y^g-i^Af^E.(-^.')'$^({Mj;ixj|T') 
= e-''^^^5:,(^=^)^<i,^({M^i(X, + di)}|r') 



^ ^ 7;'"/;" e-'-^"^'^/(^'°)(Z| -3 + 1). 



e-i-JV^E,(^=')'$^({M-iX,; + di}|T'). 



From eqn. (B16), we find that 

/(°^^"'^/2)(TZ|r) 
= (-ir)-''/2e- 



Therefore 



U{Mt)Ti = TiU{Mt). 



(B18) 



gi27r(a^i<'/3+a^K.) 



^(/3,i^-./2)(_^|_^/^) 



e 

„ I , 



= (-iT)-''/2e- 

gi27r(a'-''A'/3+a-' k) 



To obtain the algebra between U{Mt) and T2, we first 
calculate 

C/(AfT)r2$^({XJ|r) 
= U{MT)e'^^^^''"%{{X, + d2)\T) 
= e-'^^*^'(^-)'e'2.E.(^!-x.^)$^({M^iX, + d2}|r'), 

where ^2 — (0, j^). We note that 



MTd2 
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Thus we next calculate 

T2T,U{Mt)%{{X,}\t) 

= T2rie-'"^*^»(^-)'$^({Mj;ix,}|T') 

= T2e-'^^*^^(^?)'$^({M-i(X, + di)}|T') 

<f^{{M^\X,+d,+d2)}\T') 

%{{Mj;'x, + d2)}y), 

We see that 

U{Mt)T2 = e''"-T2TiL/(MT). 

Next we consider the modular transformation r — >■ 
r' = — 1/r generated by 

From 

i7(Ms)$^({X,}|T) = e-'2-^*2:,^,^x.^<i,^({M-iX,}|r'), 
we find 

C/(M5)(7(Ms)$^({X,}|t) 

= <i>^({-XJ|r). 



We see that U{Ms)U{Ms) = U{-1) generates to trans- 
formation X — —X. We can show that 

U{-~l)T^U\-l) = tI U{-l)T2UH-l) = tI (C2) 

Clearly = 1. 

Let us first calculate 

U{MsW^{{Xi\\T) 
= U{Ms)'^^{{Ms^X, + di}\T) 

Since Mgdi = ^2, we next consider 
T2U{Ms)%{{X,}\t) 

<^^{{Ms\X, + d2)}\r') 
= e-i2.iv,E.x.^^=^ci>^({M^-i(X, + d2)}|r'). 
We see that 

U{Ms)T^=T2U{Ms), 

U{Ms)T2 = tIu{Ms), 

where we have used eqn. (C2). Let us introduce T = 
e'^U(MT) and S = U{Ms), where the value of 6 is chosen 
to simplify T. We find that 

TTi = TiT, TT2 = e'''^T2TiT, 

STi = T2S, ST2 = tIs. (C3) 
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